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ABSTRACT 

The magnetic field topology plays an important role in the understanding of stellar magnetic activity. 
While it is widely accepted that the dynamo action present in low-mass partially convective stars 
(e.g., the Sun) results in predominantly toroidal magnetic flirx, the field topology in fully convective 
stars (masses below ~ 0.35 M0) is still under debate. We report here our mapping of the magnetic 
field topology of the M4 dwarf G 164-31 (or Gl 490B), which is expected to be fully convective, 
based on time series data collected from 20 hours of observations spread over 3 successive nights 
with the ESPaDOnS spectropolarimeter. Our tomographic imaging technique applied to time series 
of rotationally modulated circularly polarized profiles reveals an axisymmetric large-scale poloidal 
magnetic field on the M4 dwarf. We then apply a synthetic spectrum fitting technique for measuring 
the average magnetic flux on the star. The flux measured in G 164-31 is \Bf\ = 3.2 ± 0.4 kG, which is 
significantly greater than the average value of 0.68 kC determined from the imaging technique. The 
difference indicates that a significant fraction of the stellar magnetic energy is stored in small-scale 
structures at the surface of G 164-31. Our Ha emission light curve shows evidence for rotational 
modulation suggesting the presence of localized structure in the chromosphere of this M dwarf. The 
radius of the M4 dwarf derived from the rotational period and the projected equatorial velocity is at 
least 30% larger than that predicted from theoretical models. We argue that this discrepancy is likely 
primarily due to the young nature of G 164-31 rather than primarily due to magnetic field effects, 
indicating that age is an important factor which should be considered in the interpretation of this 
observational result. We also report here our polarimetric observations of five other M dwarfs with 
spectral types from MO to M4.5, three of them showing strong Zeeman signatures. 
Subject headings: stars: low mass — stars: magnetic fields — techniques: spectroscopic — techniques: 

polarimetric — stars: individual (G 164-31, Gl 890, LHS 473, KP Tau, Gl 896B, 

2E 4498) 



1. INTRODUCTION 

Magnetic fields in the Sun and partially convective 
stars (e.g., G, K stars and possibly early M dwarfs) 
are produced by the so-called afi dynamo mechanism 
operating at the interface between the convective enve- 
lope and the radiative core, where differential rotation is 
strongest, known as the tachocline. The action of dif- 
ferential rotation on a weak poloidal field at the base of 
the convective zone results in a large-scale and predom- 
inantly toroidal sub-surface field produced by this dy- 
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namo action (|Parkedll975f ). However, when this toroidal 
field emerges through the surface of the Sun, the mag- 
netic flux tubes are oriented in a primarially radial direc- 
tion in the photosphere (e.g. see the review by Solanki 
2009). 

Stars with masses below ^ 0.35 M0 ('^M3) are ex- 
pected to be fully convective as sugg ested by stan- 
dard models (e.g., IChabrier fc Baraffd Il997l ). This 
limit probably shifts toward low er masses due to the 
influence of the magnetic fleld (jMullan fc MacDonaldl 
I2OOII : IChabrier et all l2007| ). In fully convective stars, 
the lack of a radiative core is expected to preclude 
the aJl dynamo, maintained by the combined action 
of differential rotation {Q. effect) and cyclonic con- 
vection [a effect). This raises the question of what 
type of dynamo produces strong magnetic activity as 
observed in mid to late M an d brown dwarfs from 
photospheric Zeeman splitt ing (jJohns-Krull fc Valentil 
119961 : iReiners fc Basril I2007D to chromospheric Ha and 



Ca II IRT (iLiebert et al.l l2003t iPhan-Bao eFall 
'Schm idt et al.ll2007D . coronal X-r ay (Aud ard et al.l l 



2006 



2007 



Stelz er et all [20061: iRobrade et al. 200j), and radio 
(iBergedl2002t iPhan-Bao et all 120071: iBerger et al.l l2008t 



iHallinan et al.ll2008l: lOsten et al.ll2009f ) emissions. 

An alternative dynamo known as the so-called 
dynamo maintai ned by convection alone was proposed 
(Roberts fc Sti^ Il972l ). It has been developed for 
fully convective pre-main sequence stars (|Schiissleilll975t 
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iRudiger fc ElstoeJ [1991 iKiiker fc RiidigeJ [19971 Il999 ) 
and late M dwarfs and b r own dwarfs (iDobler et al. mm 
IChabrier fc Kiikeil 120061: iBrownind [20081). In general, 
these models predict large-scale fields, and magnetic field 
detections are now kn own for a n umber o f active objects 
in these groups (iSaanil99l iJohns-KruU fc yalentilll996l: 
Donati et all l2006al: iMorin et al.l 120081 : iBerger et all 



20091) . However, the models disagree with each other on 



the magnetic field morphology and they disagree with the 
results fro m some p hasc-rcsolvcd spc ctropola rimetric ob- 
servations ([Donati et al. 200fi_a; Mo rin et all [2008). For 
example, the latest dynamo simulations ( Browninal2008f l 
sucessfully predict that the magnetic fields possess large- 
scale fields with a lack of differential rotation in rapidly 
rotating M dwarfs. However, the simulations predict 
the field morphology in fully conv ective stars should be 
mos tly toroidal, wherea s results of lDonati et al.l (|2006al ) 
and lMorin et al.l (|2008[ ) show dominantly poloidal fields 
in 5 M3-M4.5 dwarfs. 

In this paper, we present our mapping of the M4 dwarf 
G 164-31 based on 20 hours of spcctropolarimetric obser- 
vations spread over 3 successive ni ghts. To date, there 
are only two M4 dwarfs (V374 Peg. (Ponati et al.ll2006a[ : 
G 164-31, this paper) that have been densely monitored 
over only a few rotation cycles so as to preclude intrin- 
sically long-term magnetic variability, which can affect 
the mapping result. In addition, fo r the first t i me a 
synthetic spectrum fitting techniq ue (|Marcvl[l982l : iSaail 
[19881: iJohns-Krull fc Valentil[T996h is applied for measur- 
ing the mean field strength and the filling factor in G 164- 
31. Together, the polarimetric field mapping and the 
mean field measurements provide us a better accounting 
of both large- and small-scale fields in one of these two 
rapid rotators. We also present our Zeeman signature 
search in Stokes V for 5 other M-dwarfs with spectral 
types ranging from MO to M4.5: Gl 890 (MOV), LHS 473 
(M2.5V), KP Tau (M3V), Gl 896B (M4.5V) and 2E 4498 
(M4.5V). Strong circular polarization is detected in KP 
Tau, Gl 896B and 2E 4498. 

We present our targets in § 2 and the observation and 
the data reduction in § 3. We apply the tomographic 
imaging technique for constructing the field topology of 
Gl 490B in § 4. We discuss the field topology of the 6 M 
dwarfs in § 5 and summarize the results in § 6. 

2. TARGETS 
2.1. Targets 

All the 6 targets have spectral types from MOV to 
M4.5V and their physical parameters are listed in Ta- 
ble [TJ From our mass estimates (see Table [1] and ref- 
erences therein), three M dwarfs (KP Tau, G 164-31 
and 2E 4498) have masses estimated below 0.35 Mq (or 
absolute magnitudes fainter than Mk = 6.62) and are 
expected to be fully convective according to standard 
models. Two early-M dwarfs (Gl 890 and LHS 473) are 
partially convective. The absolute magnitude of the un- 
resolved binary system Gl 896Bab is Mk = 7.28, which 
gives Mk > 7.28 for each companion. Both components 
are therefore expected to be fully convective. 

2.1.1. Gl 890 (HK Aqr), MOV 

This flari ng MO dwarf is a r apid rotator with v sim — 
70 km s~^ ([Young et al.lll98l ). A photometric period of 



0.431 days was determined bv lYoung et al.l (|1990[ ). The 
star has been well studied for magnetic activity, but no 
circu lar polarization obse rvations have been reported so 
far. iRao fc Singhl (|199(][ ) detected a steady soft X-ray 
emission from the source, indicating coronal activity. Its 
X-ray luminosity (0.2-4.5 keV) was measured by XMM - 
Newton at logLx(ergs s"!) « 29.0 (IXMM-SSQ [20Ml 
A larg e flarelike event was also detected bv ISingh et al.l 
([19991) from the X-ray light curve. Gl 890 was observed 
with the Australia Telescope Compact Array (ATCA) at 
6 cm over the full rotation phase in 1990, but the s ource 
was n ot detected at an upper fiux limit of ^0.3 mJy ([LimI 
[1991 . It was subsequently detected as a weak (below 
0.3 mJy) sl owly varyin g source with the Very Large Ar- 
ray (VLA) (|Limlll99 ^. quoting a private communication 
by M. Giidel fc A. D. Benz 1992). 

2.1.2. LHS 473 (GJ 752 A), M2.5V 

LHS 473 is an M2.5 dwarf, a companion of VB 10 (M8), 
and has a low rotati on al velocity f s im < 2.6 km s~^ 
([Delfosse et al.lll998[ ). IVaiana et all ([1981h reported a 
soft X-ray (0.2-4 .0 keV) detection with logLx = 27.11. 
iLeto et all ([20001 ) carried out VLA radio observations but 
the star was not detected in all 4 different observed radio 
bands. These observa tions indicate a low co r onal a ctiv- 
ity level in LHS 473. IStauffer fc HartmannI ([1986f ) also 
reported weak chromospheric activi ty with an Ha ab- 
sorpt io n equivalent width of 0. 35 A(|Cram fc G iampapal 
[Tip). I.Tohns-Krull fc V alentH ([l996[ ) classified LHS 473 
an inactive M dwarf. 

2.1.3. KP Tau (Wolf 1246, RX ,70339.4+2457), M3V 

This M3 dwarf, in contrast w ith LHS 473, is a rapid 
rotator with v sini = 32 km s"^ ([Mochnacki et al.l[2002 ). 
The d warf is a ROSAT source (0.1-2 4 keV ) with logLx = 
28.13 (|Fleming|ll998[ : lHunsch et al.lll999l ). implying that 
the coronal activity level in KP Tau is over an order of 
magnitude higher than that in LHS 473 (which has a 
similar spectral type). To date, no radio observations of 
the star have been reported. 

2.1.4. G 164-31 (Gl 490B), M4V 

G 164-31 is a star with a spectral ty pe of M4y 
(|Reid et al.lll995f) at a distance of 18.1 pc (|ESAII1997f) . 
which is the companion of G 164-32 (Gl 490A, M0.5V) at 
a large angular separation of 15". At the given distance 
and its apparent 2MASS if -band magnitude oiK — 8.02, 
we derive its absolute magnitude Mk — 6.73. This fast 
rotator with vsmi = 34 km s~^ exhibits high coronal ac- 
tivity in X-ray emissi on (0.3-3.5 keV) with logLx — 29.05 
([Fleming et al.l [1989( 1. No radio observations of the star 
have b een reported to date. Using the Scgra nsan et. al.l 
([2003( ) empirical relationship of Mk versus radius, we 
derive Ri, = 0.34 R©, which is i n good agr eement with 
theoretical models (e.g., Baraffc et al.lll998l) for inactive 
dwarfs (see iDemorv et al.i i200S) . This estimated radius 
yields a maximum rotational period of '^12.1 hours. 

2.1.5. Gl 896B (EQ Peg B), M4.5V 

The Gl 896 system is a quadruple Gl 89 6AabBab 
([Delfosse et al.l 119991 : [Oppenheimer et al.ll2001h . Gl 896 
Bab is a single-lined spectroscopic binary with a joint 
spectral type of M4.5 and a joint v sim of 24.2 km s~^ 
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(|Delfosse et al.l fTOQSl). Gl 89 6Bab was detected at ra- 
dio frequencies (|Giidel et aLlll993l ) of 4.9 GHz (6 cm) 
and 8.5 GHz (3.6 cm). The joint X-ray luminosity (0.1- 

2.4 ke V) of the G l 896 Aa bBab system i s logL x = 
28.84 (Giidcl ct al.] [l99l ). iRobrade et aT (2009) re- 
ported Gl 896Bab was on average a factor of 3.5 weaker 
than Gl 896Aab, we therefore derive logLx — 28.19 for 
Gl 896Bab. 

2.1.6. 2E 4498 (RX J2137. 6+0137), M4.5V 

No spectral type estimate of this dwarf is available 
in the literature. From our unpolarized (Stokes /) 
spectra obtained with ESPaDOnS (see § 3.1), using a 
spectral inde x versus spectral ty pe relationships for M 
dwarfs (e.g., iMartfn et al.lll999f i. we estimate its spec- 
tral type as M4.5ib0. 5, giving it an absolute magni- 
tude Ml K. 10.1 (e.g.. lLeggettlll992l ). Compar ison with 
the D ENIS apparent mag nitude / = 10.38 l|Epchteinl 
Il997f ) then gives a distance of ^ 11.4 pc, assuming no 
reddening. The source shows coronal activity in X- 
ray emission (0.1-2.4 keV) with a measured flux /x = 
4.52 X 10~^^ ergs s~^ cm~^ and 20 cm radio emissions 
(see iBrinkmann et al.l [2000l and references therein). At 
the source distance, we derive logLx = 28.85. The radius 
of an M4.5 dwarf is about 0.3 (jChabrier fc Baraffd 
Il997f ). resulting in a maximu m rotational period o f 

7.5 hours at v sini = 55 km s-i (jMochnacki et al.ll2002f ). 
We detected strong circular polarization in this dwarf 
(see Fig. [T|). Since its rotation period is quite short, our 
detection makes the source {Mi, ^0.2 M©) an excellent 
target for future mapping of the magnetic field on a fully 
convective star. 

3. SPECTROPOLARIMETRIC OBSERVATIONS AND 
RESULTS 

3.1. Observations 

We observed the six M dwarfs described above with 
the Canada-France-Hawaii Telescope's ES PaDOriS high - 
resolution spectrograph [R = 65,000; iDonatil [20031 ). 
which provides a wavelength coverage of 370-1000 nm. 
The spectropolarimetric mode was used to provide un- 
polarized (Stokes /) and circularly polarized (Stokes V) 
spectra. Among the six targets, G 164-31 is a particu- 
larly interesting source since it is an analog of V374 Peg 
(G 188-38), which is a rapidly rotating M4 dwarf with 
V sim » 37 km s~^ showing a large-scale poloidal field 
(|Donati et al.ll2006a^ . G 164-31 is therefore also an ex- 
cellent target on which to study the field topology in fully 
convective stars. 

Each observation consists of 4 exposures taken at dif- 
ferent polarimetcr configurations and combined together 
to filter out sp urious polarization signatures to first order 
(|Donati et al.|[l997a) . Exposure times were computed 
using the exposure time calculator for ESPaDOnS to ob- 
tain signal-to-noise ratios above 100 at the wavelengths 
of interest. This signal-to-noise value is typical of what 
is needed to d etect the circu lar polarization in M dwarfs 
(e.g., iPhan-Bao et al.ll2006f l. In the case of G 164-31, 
we monitored the star for a total of 20 hours from 2008 
March 23 to 26 with 33 observations, aiming to cover 
at least one full rotation period of the star. The full 
observing logs are given in Table [2l 

Data reduc t ion wa s performed using Libre-ESpRIT 
(|Donati et al.l Il997al) utilizing the principles of opti- 



mal extraction as given in iHornd (|1986f) . As the Zee- 
man signature is typi c ally very weak, several meth- 
ods (ISemel et al.l 120(381 iMartinez Gonzalez et all l2008t 
iDonati et all Il997af have been used in the past to in- 
crease the signal-to-noise ratio of Zeeman signatures. 
In this work, the L east-Squares Deconvolution (LSD, 
iDonati et al.l Il997a[ ) multi-line analysis procedure was 
applied to the spectra to extract the polarization sig- 
nal from a large number of photospheric atomic lines^ 
and compute a mean Zeeman signature corresponding 
to an average photospheric profile centered at 700 nm. 
Both mean Stokes / and V profiles are computed for 
all collected spectra. Figure [1] shows the LSD Stokes 
/ and V profiles of KP Tau, Gl 896B and 2E 4498 
in which we detect Zeeman signatures. The profiles of 
G 164-31 are shown in Figures [2] and [3l Using the 
basic equations listed in Wade et all (120001) (see also 
iBrown fc Landstreetlll981l : lDonati et al.lll997a[ ). we com- 
pute the net longitudinal magnetic field strength, B^, 
from our LSD Stokes / and V profiles. In the case of 
G 164-31, its time series of Stokes / and Stokes V pro- 
files are used in the following section to reconstruct the 
field topology on this M4 dwarf. 

3.2. Results 

Table [1] lists our measurements of for all targets. 
We detected strong Zeeman signatures in KP Tau, G164- 
31, Gl 896B and 2E 4498 but we did not detect them 
in the two early-M dwarfs Gl 890 and LHS 473. The 
photospheric field strength is estimated from LSD Stokes 
/ and V profiles, and the chromospheric field is computed 
from the Stokes profiles of only the Ha line. 

For G 164-31, we find that the Zeeman signatures re- 
peat after 0.540±0.006 day (12.96±0.14 hours), which is 
in the expected range of its rotation period. We thus con- 
clude that the rotation period of this dwarf is 12.96 hours. 
We note that this value is slightly greater than the es- 
timated maximum rotation period of 12.1 hours for this 
star. The difference is possibly due to G 164-31 hav- 
ing a larger radius than the expected value from the- 
ory. Matching unpolarized (Stokes /) spectra of G 164- 
31 with that of slowly rotating references requires a pro- 
jected rotation velocity of tisim = 41 zfc 1 km s ~^. Our 
value is greater than one of lFleming et al.l (|1989f ). vsini — 
34 ± 8.5 km s^^, though both measurements are consis- 
tent to within the uncertainties. Dues to its smaller un- 
certainty, we adopt our vsini = 41 km s~^ for G 164-31. 
/^From Prot and vsini, we derive Rsini = 0.44 ± 0.02 R©. 
Our i?sini is ~ 30% larger than the estimated stellar ra- 
dius (see § 2.1.4). It is unlikely that i exactly equals 90°, 
and the results obtained in the imaging process decribed 
below are similar for 60° < z < 80°. We thus set z = 70°, 
which results in the G 164-31 radius being '^38% larger 
than the typical value for M dwarfs of the same lumi- 
nosity. This is likely due to the fact that G 164-31 is a 

^ We used a line list l IKurucj 119931 ) corresponding to M spec- 
tral types matching that of the six dwarfs. Approximately 5000 
intermediate to strong atomic spectral lines with an average Lande 
factor g^ff of ~ 1.2 are used simultaneously to retrieve the average 
polarization information with typical noise levels of RiO.08% (rela- 
tive to the unpolarized continuum level) per 1.8 km velocity 
bin and per individual polarization spectrum. This represents a 
multiplex gain in S/N of about 10 with respect to a single line 
analysis. 
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young object belo nging to a moving cluster identified by 
lOrlov et al.l (|1995l ). We refer the reader to § 5 for furtlicr 
discussion. 

We monitored G 164-31 for 20 hours spread over three 
successive nights, aiming for coverage over a full rota- 
tional cycle. The true period (12.96 hours) is a bit longer 
than our expected maximum value (12.1 hours). With 
this period, our observations cover only half of the ro- 
tational phase. In other words, only half of the stellar 
surface was monitored. A few of the Stokes / and V LSD 
profiles are discarded due to the contamination of moon- 
light or flaring events. Figure [2] and [3] show Stokes / and 
V LSD profiles of G 164-31. 

4. MODELING OF MAGNETIC FIELD TOPOLOGY OF G 

164-31 

To reconstruct the map of cool spots and magnetic 
fields at the surface of G 164-31 fro m the set of observed 
LSD Stokes / and V, we use the iDonati et ahl ()200lD 
magnetic mapping code which employ es the Zeeman - 
Doppler imaging (ZDI) technique (Sc mel et all Il989f ). 
The typical longitude resolution achieved at the equa- 
tor is « 15° or w 0.04 rotation cycle. The basic idea 
of the code has be e n fully described in s evera l papers 
(iBrown et all fl99H IDonati et al.l Il997bl [200lL l2006bl : 
iMorin et al.ll2008( r^In this paper, we briefly describe the 
principles of the modeling work and refer the reader to 
the above papers for further details. 

4.1. Mapping spots 

Cool spots on the surface of a rapidly rotating star pro- 
duce distortions in the stellar spectral lines. If the star is 
rotatin g fast enough (at least 20-30 km s~^ for late- type 
dwarfs. TVogt et al.lil987| ). the shape of the star's spectral 
line proflles are dominated by rotational Doppler broad- 
ening. In this case, there is a strong correlation between 
the position of any distortion within a line profile and the 
position of the corresponding spot on the stellar surface. 
A spectrum of the rapid rotator is a one-dimensional im- 
age that is resolved in the direction perpendicular to the 
stellar rotation axis and the light-of-sight. By observing 
the star at different rotation phases, a two-dimensional 
image of the spot distribution can be reconstructed us- 
ing the so-called Doppler imaging (DI) technique (see 
IVogt et all Il987l and references therein). In the imag- 
ing process, the stellar surface is divided into a grid of 
10 00 elemen t ary ce lls. Using the spot-occupancy model 
of ICameronI (|1992f ) , the local line profile at each grid 
point of the surface is described as a linear combination 
of two reference proflles, one representing the quiet pho- 
tosphere and one for cool spots. Both reference profiles 
are assumed to be equal and only differ by their relative 
continuum levels. For the assumed reference profile, we 
can use either the LSD profile of the very slowly rotating 
inactive M4 dwarf Gl 402 or a simple Gaussian profile 
with similar full width at half-maximum and equivalent 
width. Both options yield very similar results, indicat- 
ing that the exact shape of the assumed local profiles has 
very minor impact on the reconstructed images, as long 
as the rotational velocity of the star is much larger than 
the local profile. Each element on the surface is quan- 
tified by the local fraction occupied by spots. The spot 
occupancy ranges from (no spots) to 1 (complete spot 
coverage). The maximum entropy image reconstruction 



technique (jSkilling fc Brvanlll984l : IVogtlll980[ ) is used to 
find the image having the smoothest variation (or in some 
sense the simplest image), giving the least contrast be- 
tween spots and the quiet photosphere. 

4.2. Mapping magnetic fields 

The field is described as the sum of a poloidal and 
a toroidal compone nt, both expressed as spherical har- 
monics expansions (iJardine et al.lll999( ). For a given set 
of the complex coefficients of the spherical harmonic ex- 
pansions, one can produce a corresponding topology (or 
a Stokes V data set) at the stellar surface. The code 
uses the maximum entropy image reconstruction tech- 
nique, in which entropy (i.e., quantifying the amount of 
reconstructed information) is calculated from the coeffi- 
cients of the spherical harmonics expansions. The imag- 
ing process starts from a null magnetic field and itera- 
tively adjusts the spherical harmonic coefficients in order 
to match the synthetic Stokes V profile to the observed 
LSD one. The code uses a multidirectional search in the 
image space until the required maximum entropy image 
is obtained. This corresponds to an optimal field topol- 
ogy that reproduces the data at a given level, i.e., 
usually down to noise level. Since the inversion problem 
is partially ill posed, we use the entropy function to select 
the magnetic field map with lowest information content 
among all those reproducing the data equally well. 

To calculate the synthetic Stokes V profiles for a given 
field topology, the surface was again divided into a grid 
of 1000 elementary surface cells. In each cell, the three 
components (radial, azimuthal and meridional in spheri- 
cal coordinates) of the vector field are estimated directly 
from the spherical har monic expansion. The code uses 
the a nalytic solutions (jLandolfi &: Landi Degl'Innocentil 
|1982|) of Unno-Rachkovsky's radiative transfer equations 
to calculate the contribution to the Stokes V profiles of 
all visible cells at each observed rotation phase. The free 
parameters in the Unno-Rachkovsky equations are ob- 
tained by fitting the LSD Stokes I profile of a very slowly 
rotating and weakly active star with a similar spectral 
type (e.g., Gl 402). To obtain the best fit (in both am- 
plitude and width) between the synthetic and observed 
Stokes V profiles, the code introduces a filling factor fc, 
which represents the fractional amount of circular fiux 
being constant over the whole stellar surface, and specu- 
lates that large-s cale fi elds in M dwarfs can be structured 
on a small scale (|Donati et al.ll2008l ). The circularly po- 
larized fiux from each cell is thus /cl4oc, where Vioc is 
the Stokes V profile derived from the analytic solutions 
of the Unno-Rachkovsky equations. 

Differential rotation is also implemented in the calcula- 
tion process of the synthetic Stokes V profiles for yield- 
ing the best fits to the observation. For this purpose, 
we assume the rotation rate varies with latitude, 6, fol- 
lowing a solar-like differential rotation law as r2(0) = 
^^cq — sin^0, where ilcq is the angular rotation rate at 
the equator, and dfl is the difference in angular rotation 
rate between the equator and the pole. Differential rota- 
tion is detected when of the fit to the data shows a well 
defined minimum in the range of f2cq and dQ values. In 
the case of G 164-31, we did not obtain a clear minimum 
in the explored Qcq — dfl range. This indicates that our 
data are not suitable for measuring differential rotation 
given the fairly simple (and low amplitude) rotational 
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modulation of the Stokes V profiles we observe. This is 
mostly due to the fact that the field distribution lacks 
well defined features at different latitudes from which 
differential rotation can be estimated. Bas ed on the pre- 
vious observations of analogs of G 164-31 ( Donati et al.l 

[2b06a: .Morin et al 2008 1 , we therefore assumed that this 

star rotates as a rigid body while performing the field 
mapping. The obtained map of G 164-31 is shown in 
Figure [H 

4.3. Results 

Figure |4] presents our maps of the spot occupancy and 
the magnetic field at the stellar surface of G 164-31. For 
the half stellar surface that was not observed, the field 
topology has been derived using a spherical harmonic 
expansion based on the collected Stokes V data set. Our 
brightness map shows a low contrast spot at the pole 
of the star. The magnetic field morphology on G 164- 
31 shows that the radial field component is dominant, 
which appears similar to what is seen in V374 Peg. The 
azimuthal and meridional components are negligible. 

The average large-scale flux is 680 G and the poloidal 
field dominates with 99% of the energy content. Most 
(95%) of this poloidal field is stored in the mode cor- 
responding to a dipole aligned with the rotation axis 
(spherical harmonics degree I = 1 and order m = 0). 
One should note that due to the high vsinz of G 164-31, 
for any value of the filling factor < /c < 1, wc find 
similar results. Hence, we arbitrarily set /c = 1. We 
also examine the binarity of G 164-31. Based on the 
LSD Stokes / profiles, we find an average radial velocity 
value of —6.7 ± 0.2 km s~^. No significant variation in 
radial velocity is found. Therefore, our current data does 
not reveal any hint of an additional close-in component 
around G 164-31. 

5. DISCUSSION 

In this section, we first briefly discuss the magnetic 
field properties of the stars whose fields are not mapped. 
These stars are Gl 890, LHS 473, KP Tau, Gl 896B and 
2E 4498. We then focus the discussion on the field topol- 
ogy in the photosphere and the chromosphere of G 164- 
31. 

5.1. Individual stars: Gl 890, LHS 473, KP Tau, Gl 
896B and 2E 4498 

We did not detect Zeeman signatures in Gl 890 and 
LHS 473. For the latter case, the non-detection in both 
the photosphere and the chromosphere indicates that 
LHS 473 is inactive, likely due to its low wsinz, and this 
result is consistent with previous observations showing 
its low coronal activity level. However in the fastly ro- 
tating case of Gl 890, the non-detection suggests that 
the field topology on this star may be different from 
what has been observ ed in more slowly-rotating early-M 
dwarfs (Donati et al, 2008), whose fields are large-scale 
and dominantly toroidal producing Stokes V signatures 
detectable at any time. Since the detections of X-ray and 
radio emission indicate that Gl 890 is magnetically ac- 
tive, phase-resolved spectropolarimetric observations of 
this star are therefore needed to clarify its magnetic field 
topology 

For the three remaining M dwarfs, all these fast rota- 
tors show Zeeman signatures, and they are particularly 



strong in Gl 896B and 2E 4498 (Fig. [T]). This result in- 
dicates the presence of strong, large-scale fields on these 
stars. Since KP Tau (0.31 Mq) and 2E 4498 (0.21 Mq) 
are fully convective, they are good targets for studying 
the field morphology by carrying out spectropolarimet- 
ric a nd simultaneous rn ultiple-wavelength observations 
(e.g., iBerger et aLll2009D . We note that Gl 896Bab is 
an unresolved binary. Therefore, its Zeeman signature is 
contributed by two components if both are magnetically 
active. 

5.2. G 164-31 

The G 164-31 magnetic field topology, which is mostly 
poloidal with most of the magnetic energy concentrated 
within the lowest order axisymmetric mode s, is very sim- 
ilar to what has been observed in V374 Peg (lEjonati et all 
I2006a[ ) and other mid-M dwarfs (jMorin et al.ll2008l ). The 
average large-scale magnetic flux on G 164-31 is 0.68 kG, 
which is significantly smaller than the previous measure- 
ments of few kiloGauss magnetic fluxes in active mid-M 
dwarfs using synthetic spectrum fitting techniques (e.g., 
iJohns-KruU fc Valentil lT996') . This is possibly due to the 
ZDI technique being sensitive to only large-scale and 
simple structures, while Zeeman signatures from mag- 
netic regions with complicated or small-scale topology 
may cancel each other in circularly polarized spectra 
whereas these signatures add up in unpolarized spec- 
tra. To explore this "missing" magnetic flux, we used 
the synthetic spectrurn fitt ing technique described in 
iJohns-KriiU fc Valentil (fl996h to measure the mean mag- 
netic flux on the surface of G 164-31. Briefly, spectra 
in the wavelength interval around the Zeeman-sensitive 
Fe I line at 8468.40 A are analyzed. In stars as cool as 
G 164-31, this line is actually blended with a Ti I line at 
8468.47 A which is similar in strength and also Zeeman 
sensitive. The spectrum synthesis models both lines plus 
numerous TiO lines in this wavelength region. Because 
of the ubiquitous TiO, spectral changes due to magnetic 
fields are seen more clearly in the ratio of active to in- 
active line profiles. Therefore, spectra of G 164-31 were 
divided by inactive references. In our case, we used two 
inactive reference stars GJ 725B (or LHS 59, M3.5V) and 
GJ 876 (or LHS 530, M4V) which are meant to be iden- 
tical in all respects (T^ff, logg, [M/H], etc.) to G 164-31 
except for the field. Line profiles were synthesized with 
and without magnetic fields (see Johns-KruU & Valent^ 
Il996l and references therein) and the ratio of active to 
inactive profiles was fit by adjusting the magnetic field 
strength, B, and its filling factor, /, in the model spec- 
trum of G 164-31. The line profiles are generated as- 
suming a uniform radial field everywhere with a uniform 
filling factor. 

The best fits were obtained with Bf — 3.1 kG and 
Bf = 3.3 kG for GJ 725B and GJ 876. The fit to the 
observations using GJ 876 as the inactive reference star 
is shown in Figure 5. Due to the high vsmi of G 164-31, 
we can not measure B and / separately. The model fit in 
the line core is not as good as we would like, and is only 
slightly improved using GJ 725B as the reference star. 
We note though that if the radius of G 164-31 is indeed 
larger than expected due to youth (see below), its grav- 
ity is likely somewhat lower than that of cither inactive 
comparison star. Our spectrum synthesis indicates that 
the 8468 A feature should actually weaken somewhat at 
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lower gravity, which would result in the discrepancy in 
the core getting worse by ^ 0.005 in the ratio. A possible 
way to correct for this would be to fit a distribution of 
magnetic field strengths on the stellar surface as has been 
done for M dwarfs (Johns-KruU & Valenti 2000) and on 
T Tauri stars (e.g. Johns-KruU et al. 1999). The current 
fit (Fig. 5) attempts to use a single field strength to fit 
both the wings and the core, whereas a distribution of 
magnetic fields gives more flexibility in fitting the entire 
profile and can result in bigger changes in the line equiv- 
alent width than produced when using only a single field 
value. We therefore adopt Bf = 3.2±0.4 kG for G 164-31 
as the best estimate from our single field fitting described 
above, noting that this may be a slight underestimate. 
This value for Bf is larger than the one measured from 
the ZDI technique by a factor of 4.7. Our result is consis- 
tent with the previous measurements from t he literature 
(e.g. EV Lac, iJohns-Krull fc Valentil ITOOl iMorin et all 
|2008| ). suggesting that in active mid-M dwarfs a signifi- 
cant part of the magneti c energy is stored insmall-scale 
structures. As found bv iReiners fc Basrl (|2007f) . it ap- 
pears that a large majority of the magnetic energy is 
stored in the small scale field of G 164-31. Apparently, 
a successful dynamo model should be able to produce 
both large-scale poloidal fields and small-scale features 
in rapidly rotating mid-M dwarfs. 

To study chromospheric activity, we computed the Ha 
emission equivalent widths for each exposure. Our Ha 
emission light curve is shown in Fig. [6l One should 
keep in mind that since each observation to obtain a 
Stokes V profile consists of 4 exposures, measuring Ha 
emission at individual exposures therefore increases time 
resolution of the light curve. Figure [6] also indicates 
that the Ha equivalent width variations are likely sinu- 
soidal and they are modulated by the stellar rotation 
with Prot = 0.54 day determined from the Stokes V 
profiles. This effect h as also been observed in late-M 
and brown dwarfs fsee lBerger et al.|[200& and references 
therein). Flaring events were seen strongly in the last 
sequences of the first observing night and moderately in 
the middle sequences of the second night. These events 
were also observed in the H/? emission profiles. While 
the likely rotational modulation and nonzero minimum 
of the Ha light curve imply either a large-scale (poloidal 
or toroidal) chromospheric field or a localized concen- 
tration of small-scale magnetic activity, it is reasonable 
to assume that the chromospheric field in G 164-31 has 
both a large scale poloidal component and small-scale 
structures as observed in the photospheric field. This 
field configuration probably dominates in both the pho- 
tosphere and the chromosphere of the star. We note that 
simultaneous observations at multiple wavelengths of the 
source will place more contraints on the chromosphe ric 
and coronal field morphology (e.g.. iBerger et al.]l2009| ). 

It is very interesting to note that the radius of G 164- 
31 is at least 30% larger than that estimated from 
standard models for M dwarfs at the same luminosity 
(see § 2.1.4). There are two possible interpretations 
of this discrepancy. First, the magnetic field effects of 
reduced convective efficiency due to fast rotation and 
lar ge field strength s , and/ or to spot coverage as proposed 
bv IChabrier et all (12007 *) might yield a cooler effective 
temperature and thus a larger radius in low-mass 
stars, keeping their luminosity unchanged {L oc T^^B?). 



Our spot mapping (Fig. [4]) shows the fraction of the 
stellar surface covered by the low contrast spot is rel- 
atively small in G 164-31. However, one should note 
that the imaging technique is only sensitive to spots or 
spot groups with sizes comparable to our resolution el- 
ement. We therefore can not rule out the possibility of 
unresolved small spots spread everywhere on the stellar 
surface. Hence, we can not conclude whether the effect 
due to spot coverage is s ignificant in this star or not . Ob- 
servational results (see 'RibasI 120061: [Morales et al]|2008l 
and references therein) suggest that these effects might 
yield a 10 — 15% larger radius in G 164-31, there- 
fore the discrepancy in radius of over 30% can not be 
explained by only those effects. Second, we explore the 
possibility that G 164-31 is a young M4 dwarf, its radius 
thus larger than one computed for old M dwarfs. We 
note that the if-band absolute magnitude of G 164-31, 
Mk = 6.73, is 0.67 mag brighter than the average value 
for M4 dwarfs"'^'^, supporting the young nature of the star. 
From a literature search, the binary system G 164-31 (Gl 
490B)-|-G 164-32 (Gl 490A) ind eed belong s to a moving 
cluster previously identified by lOrlov et a l. ( 199E j) (for 
a rev iew on young nearby stars, see IZuckerman fc SongI 
I2004f) . One should note that the primary component 
G 164-32 also exhibits high coronal X-ray emission with 
log(Lx/iboi)= —3.23 (Fleming et al. 1989.) and ch romo- 
spheric Ha emission (|Stauffer"^Hartmannl 1 1981 ) . indi- 
cating activity typical in young early-M dwarfs. We mea- 
sured an upper limit for the Li A6708 equivalent width of 
17 mA for G 164-31. This gives the star an age older than 
10 Myr since early-M dwarfs (>M4) are expected to com- 
pletely deplete their Lithium within ^10 Myr. We sug- 
gest G 164- 31 is an ana log of the M4 dwarf HIP 112312A 
at 23.6 pc (|ESAl[l99l . which is a member of the ~12 
Myr old (3 Pictoris moving group (jSong et al.ll2002[ ). To 
estimate the m ass and age of G 16 4 -31, w e used the 
IChabrier et"aI1 (|2007f l: iBaraffe et all (|1998[ ) theoretical 
models with constraints of Mk = 6.73, I — J = 1.61, 
and the deduced radius i?* — 0.47 R© of G 164-31, and 
we also considered the magnetic field effects that might 
yield about a 10-15% larger radius. We then derived 
Mi, ^ 0.15 Mq at an age of about 25-30 Myr, mak- 
ing G 164-31 the least massive M dwarf whose magnetic 
field has been mapped to date. More observations are 
needed to precisely determine the fundamental parame- 
ters of this moving cluster. 

6. SUMMARY 

Based on our circularly polarized and unpolarized spec- 
troscopic observations, using both tomographic imaging 
and synthetic spectrum fitting, we reveal the mainly ax- 
isymmetric large-scale poloidal magnetic field and the 
small-scale field structures storing a significant portion 
of the magnetic energy in the photosphere of the M4 
dwarf G 164-31. The modulation of the Ha emission 
light curve suggests that the field in the chromosphere is 
stable and possibly poloidal like that in the photosphere. 
Our detection of circular polarization in the single and 
rapidingly rotating M dwarfs KP Tau (M3V), 2E 4498 

We sele c ted 9 single M4 dwarfs within 10 pc listed in 
IDelfosse et al.l l|1998l ) and computed their K-hanA absolute mag- 
nitude from trigonometric parallaxes and 2MASS /f-band mag- 
nitudes available in the VIZIER, database. This work yielded an 
average value Mk = 7.4. 
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(M4.5V) make them good targets for mapping of the field 
morphology in mid-M dwarfs in the future. 
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TABLE 1 

Physical parameters of the 6 M dwarfs 



Star 


SpT'' 




logLboi'' 




log(i?x)° 






Ha EW 


Bz(Hq) 








{ergs s~-^) 


(km s~-^) 






(G) 


(A) 


(G) 


Gl 890 


MO 


5.41 


32.26 


70 


-3.26 


0.57 


<18 


0.64±0.01 


<22 


LHS 473 


M2.5 


5.83 


32.09 


<2.6 


-4.98 


0.48 


16±7 


-0.29±0.02 


<30 


KP Tau 


MS 


6.82 


31.66 


32 


-3.53 


0.31 


-104±32 


0.73±0.05 


-580±220 
















-92±32 


0.73±0.05 


210±110 


G 164-31 


M4 


6.73 


31.66 


41 


-2.61 


0.33 


680S 


3.15±0.0is 


530±60S 


Gl 896B 


M4.5 


7.28 


31.45 


24.2 


-3.26 


0.25 


296±40 


3.5±0.3 


-150±60 
















294±37 


3.6±0.3 


110±60 
















279±38 


4.8±0.5 


160±15 


2E 4498 


M4.5 


7.60 


31.31 


55 


-2.46 


0.21 


-440±84 


4.7±0.5 


-270±100 



Spectral type from the literature. 

Absolute ii'-band magnitudes computed from trigonometric parallaxes (Gl 890, LHS 473, G 164-31, Gl 896B) or 
spectroscopic distances (KP Tau, 2E 4498) and 2MASS A'-band magnitudes available in the VIZIER database. 

B olometric lumin osity estimated using the bolometric correction BCk versus i^-band magnitude relationship in 
Prinnev et all ||T993) . 

Rotational velocity, see text (§ 2.1) for references, except G 164-31 (this paper). 

" -Rx = ix/^bol- 

' Stellar masses estimated using the Mk versus mass empirical relationship in lDelfosse et all (120001 ). 
^ For G 164-31: the value listed for an average magnetic flux Bfc computed using the ZDI technique (see § 4.3), 
assuming /c = 1. Its H« equivalent width and chromospheric field strength B^, are estimated at cycle E = 0.9181 
where the Zeeman signature appears very strong (see Fig.[3ll. 
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TABLE 2 
Observing logs for the 6 M dwarfs 



Target 


Date 


HJD^ 


UT 


Exposure 


Peak of 


Cycle 






(2453000+) 


(h:m:s) 


time (s) 


S/N 




Gl 890 


9007 '^on 9Q 


1 '^79 Q1 


nQ-4Q-1 


4x200 


166 






900^1 Son 1^^ 


7-d70Q 




4x300 


511 




KP Tau 


9007 Son 9Q 


-LO / O . UO -LO / 


1 '^■nS-9zL 
j-O . uo . 


4x500 
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^ Heliocentric Julian date (UTC). 

^ Rotational cycles E are computed using ephemeris HJD = 2454549.54 + 0.54E, 
only available for G 164-31. 
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Fig. 2. — LSD Stokes I profiles of G 164-31 (black line) along with the maximum entropy fits (red line) to the data. The rotational 
phases are indicated to the right of each profile. 
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Fig. 3. — LSD Stokes V profiles of G 164-31 (black line) along with the maximum entropy fits (red line) to the data. The rotational 
phase (right) and the 3cr error bar (left) of each profile are shown accordingly. 
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Fig. 4. — Magnetic field topologies of G 164-31, reconstructed from the series of LSD Stokes V profiles. The radial field {bottom) 
dominates in the dwarf (intensity scale in Gauss), the azimuthal and meridional component are negligible. The spottedness (top) at the 
stellar surface is also reconstructed, the color scale represents spot occupancy for the dwarf (l=complete spot coverage). The dwarf is 
shown in flattened polar projection extending down to latitudes of —30°. The equator is described as a bold circle. Radial ticks outside 
the plots indicate the phases at which the dwarf was observed, a half of the dwarf surface had properly been monitored. 



Magnetic Fields in M Dwarfs 



15 



1.04 



100 

00 

O 



tn 

o 



0.92 



— I 1 1 r- 



— I 1 1 1— 



— I 1 1 r- 



8460 



_i I I i_ 



_l I I L_ 



8465 8470 ^ 

Wavelength (A) 



_i I I i_ 



8475 




Fig. 5. — Best-fit modeling ratio of active to inactive profiles (thick line). The observed line profile ratio constructed by dividing the 
spectrum of G 164-31 (Gl 490B) in the wavelength region around the magnetically sensitive Fe I line at 8468.4 A by that of the inactive 
M4 dwarf Gl 876 is shown (thin line). 
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Fig. 6. — Ha equivalent width light curve. The curve appears likely sinusoidal, a sinusoidal curve (dashed-line) is also shown for 
a comparison. The first vertical tick at 0.19 in the HJD oflfset (HJD— 2454549.79717) position indicates the bottom of the light curve 
obtained in the first night. The remaining ticks are drawn with a distance of 2 X Prot between two consecutive ticks where the star is 
observable, where I\ot = 0.54 days is the rotational period determined from the Stokes V profile set (see § 3.2). The observed bottoms of 
the light curve likely repeat after 2Prot , suggesting that the H« emission is modulated by the stellar rotation. Flaring events were observed 
in the first (strong) and second night (moderate). The decay and enhancement of Ha emission were also observed in the third night at 
around HJD offset=2.1 and 2.35, respectively. This is possibly due to a decrease (or an increase) of the number of active regions. Changes 
in the Ha (also H^S) line profile shape (e.g., smoothness) seen in our spectra support this possibility. 



